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Abstract 

One problem of wide interest involves estimating expected crossing-times. Several tools have been 
developed to solve this problem beginning with the works of Wald and the theory of sequential anal- 
ysis. An extension of his approach is provided by the optional sampling theorem in conjunction with 
martingale inequalities. Deriving the explicit close form solution for the expected crossing times may be 
difficult. In this paper, we provide a framework that can be used to estimate expected crossing times 
of arbitrary stochastic processes. Our key assumption is the knowledge of the average behavior of the 
supremum of the process. Our results include a universal sharp lower bound on the expected crossing 
times. 



Let {A t } t >o be a non-negative, measurable process. Set q = sup < s < t 
t > and T r = inf{£ > : Xt > r}, r > 0. In this paper, we present bounds on the expected hitting time 
of Xt via a nonrandom function a(t) := E [<;(£)] with a~ (s) = oo if a(oo) < s. In particular, we derive 
the following sharp universal lower bound , E[a(T r )] > r/2 and if in addition a(t) is concave, we have 
■Epr] — a_1 ( r /2)- This bound has the optimality property that in the case of non-random continuous 
processes, the following identity is valid: T r — a (r). Therefore, (up to the constants) the method 
provides a unified approach to boundary crossing of non-random functions and stochastic processes. 
Furthermore, for a wide class of time-homogenuous, Markov processes, including Bessel processes, we 

are able to derive an upper bound E[a(T r )] < 2r, which implies that sup r 

inequality motivates our claim that a(t) can be viewed as a natural clock for all such processes. The 
cases of multidimensional processes, non-symmetric and random boundaries are handled as well. We also 
present applications of these bounds on renewal processes in Example 10 and other stochastic processes. 
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1 Introduction 

One problem of wide interest in the study of stochastic processes involves estimating E[T r ], the expected 
time at which a process crosses a boundary. Following the work of [7], we let X t , t > 0, be a measurable 
process with the first passage time T r = mf{t > : X t > r}, r > at the level r. Deriving the explicit closed 
form solution for E[T r ] can sometimes be difficult and so we arc interested in finding the possibility of ob- 
taining bounds for E[T r ] above or below by some functions that are related to a(t) := -E[sup 0<;j<t X s ]. This 
approach was introduced in de la Pena [7j (published in Section 2.7 of [10]) in which the author constructed 



bounds for E[T r ] for a class of Banach- valued processes with independent increments via decoupling. The 
bounds derived are of interest in applications where the moments of the maximal process can be readily 
obtained. 



The main idea consists of an extension of the concept of boundary crossing by non-random functions to 
the case of random processes. a(t) can be intuitively interpreted as a natural clock for all processes with 
the same a(t). Here, we assume that we have information on a(t). Furthermore, we assume that X t is a 
measurable separable process, sup 0<s<t X s is, therefore, well-defined and T r is measurable. 

We would like to draw an attention to readers that the results derived in this paper provide a decoupling 
rcintcrpretation (in the context of boundary crossing) of the results due to Wald that concern randomly 
stopped sums of independent random variables (see PS]), to Burkholder and Gundy on randomly stopped 
processes with independent increments (see [5]) and to Klass on bounds for randomly stopped Banach space- 
valued random sums (see [13] and [M]). 

Recall that for {Xi} a sequence of iid random variables adapted to Ti = <x(Xi, ■ • ■ , Xi) with i?[Xi] = fi, 
where < oo, we have 

E[S T ] = fiE[T], 

where T is a stopping time adapted to a{Xi}, S n — Y^i=i an d E[T] < oo. Furthermore, if £7[-Xi] = 
and E[Xf] < oo, then 

E[S T } = E[Xf]E[T], 

whenever E[T] < oo. These celebrated results are known as Wald's first and second equations. As we can 
observe, St = S£i XiI(T > i) is composed of summands XiI(T > i) that are products of independent 
variables. This motivate the idea of replacing {XiI(T > i)} by {XiI(T > i)}, where {Xi} is an independent 
copy of {Xi} and independent of T as well. As a result, the above Wald's equations can be viewed from the 
decopuling perspective where, if we denote S n = J2i=i Xi, then, whenever E[T] < oo and -E[AT^] < oo, 

E[S T ] = E[S T ] = E[S f ] and E[S T ] = E[S T ] = E[S 2 f ]. 

It is important to realize that the variables St and St can have drastically different behavior. One may 
consider a case in which Xj's are iid random variables with Pr{Xi = 1} = Pr{X; = —1} = |. Define 
T = inf{n > : S n = a or S n — —b} for some integers a, b > 0. It is easy to see that St can only take a 
value of either a or —b whereas Sy's value is not restricted to these two choices. 



An extension of Wald's second equation to the case of independent random variables having finite second 
moments with T as a stopping time defined on the X's was studied by de la Peha and Govindarajulu |11) . 
The following bound is sharp: 

E[S T ] < 2E[S T ], (1) 

where S n = Yh=i Xi and with S an independent copy of S. If we let T r = inf {n : S 2 > r} and a(n) = 
£[maxKj< n Sj], then 



max S ? 

l<j<T r J 



r < E[S Tr ] < 2E[S Tr ] < 2E 
hence 

E[a(T r )] > r -. 

which is closely related to our main result; see Proposition 1 below. 



Along this vain, Klass [13] obtained results on a best possible improvement of Wald's equation. In 
his work, Klass obtained bounds for stopped partial sum of independent random elements taking values if 
Banach space (B, \\ ■ ||). To be specific, he derived that 



E 



max $(\\S n \ 

Kn<T 



< 20(18 a )£; 



max $(\\S n \ 

Kn<T 



2 



where <&(•) : [0, oo) is a nondecreasing function such that $(0) = and for a > 0, <&(ca;) < c Q $(x) for all 
c > 2, x > . Furthermore, in |14j . the corresponding lower bound was also derived and hence 



c a E 



max $(||SU|) 

Kn<T V " " y 



< £ 



max $(||SJ|) 

Kn<T V " 



< 20(18 a ).E 



max 

Kn<T V " ' 



(2) 



The corresponding counterpart for processes defined in continuous time domain was developed by [9] (see 
also [5]) in which they showed 



c n E 



sup $(\\N t \\) 

Kt<T 



< E 



sup $(||JV t ||) 

Kt<T 



< 20(18")S 



sup m\N t \\) 

Kt<T 



where Nt, t > is a B- valued process continuous on the right with limits from the left with independent 
increments with (B, || ■ ||) is a separable Banach space. As observed above, the inclusion of a(t) facilitates 
the decopuling between the random stopping time T and the underlying process X t . All these results are 
closely connected with sequential analysis, for details, readers are referred to Lai [15] for a recent survey. 

In our results, below, we use E[a(T r )] as a key quantity. Here, we review our interpretation of this 
quantity. If we define 

M(t) — sup A(s), and X t is continuous 

0<s<t 

then, M(T r ) = r with probability one. But, in contrast, 

E[a(T r )} := [ {E[M(t)])dF Tr (t) = [ (E[M(t)])dF Tr (t), 



which can be interpreted as E[a(T r )] = E[a{T r )\ = E[M(T r )], where T r T r with T r and T r independent. 
Thus it fits into the decoupling theme discussed above. This is further discussed in the remarks following 
Proposition 1. 

In this paper, we devlope upon [4] and [8] to provide a universal sharp bound for E[a(T r )} > |, which 
under a concavity assumption on a(t) gives E[T r ] > a _1 (|); compare de la Pena and Yang 
following bound is presented: 

E[T r ] > (1 - e)a- 1 (er), e 6(0,1). 
In addition, for a wide class of processes, we show that 



in which the 



(3) 



l U<E[T r ]<a-\2r) 



E[a{T r )} 



as well as the stability property sup 

r>0 

without the concavitiy assumption on a, gives: 

1 

2° ^2 



< 1. The above result, coupled with Eq. (1.1) of [12], 



(-)<E[T r ]<a-'(2r). 



It should be noted that, the theory of first passage times for random processes has been extensively 
developed in recent times. In particular, the distribution of the first hitting times of Brownian motion has 
been studied through inverse Gaussian distribution; see [19]. A similar approach to the first hitting times 
involving Levy processes is also available. Readers may refer to [T5] for details. 

The typical methods, however, assume full knowledge of the distribution. In contrast, our approach pro- 
vides bounds for all processes with a common a(t) = .E[sup 0<s<t X s ] based on the approximate knowledge of 
moments of the maximal process, or i?[A t ]. Even in situations when the distribution of the process is known, 
the quantity E[T r ] might not be easily obtained as shown in Example 7 in which we study the relative growth 
of the boundary crossing of a three-dimensional Brownian motion and related processes. (Renewal processes) 

The rest of the paper is organized as follows: In section 2, we obtain upper and lower bounds on E[T r ], 
as well as bounds on E[a(T r )]. Section 3 elaborates some possible extensions of our methodology that can 
handle siutations in which the expected first hitting time is hard to obtain. An application of the bounds 
derived are presented in Section 4, followed by Section 5, which summarizes the paper. 
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2 Main results 



With the above definitions, let a — inf{f > : a[t) > £} and a = oo if a(t) < £,Vt. We have the 
following proposition. 



Proposition 2.1 For all non-negative, measurable process X t with 
E [sup < s < t X s ] = a(t), 

\ < E[a(T r )} 

and the bound is sharp. Furthermore, if a(t) is assumed to be concave, we obtain 

a-\ r -)<E[T r ]. 

A decoupling reinterpretation of Q is given as follows: Assuming that X t is continuous, we have 



(4) 



(5) 



E 



sup X s 

0<s<T r 



E[M(T r )}, 



we can rewrite Q as 



1 



-E 









sup X s 


< E 


sup X s 


0<s<T r 




0<s<T r 



2 2 

in which case the underlying process and its random stopping time are independent (decopuled). 
Proof 1 Let F be the cdf of T r . If F is continuous, then 



sup X(s) 

0<s<t 



a(t) = E 



> rPr< sup X(s) > r 
[a<s<t 

= rPr{T r <t} = rF(t). 

Due to the continuity of F, it is easy to see that F(T r ) - U(0, 1) and E[F Tr {T r )} = |. As a result, we can 
conclude that E[a(T r )] > |. More generally if the distribution of T r is not necessarily continuous then since 



E[F(T r )} = Pr (T r < T r ) 



= Pr ( T, = T, 
1 



1 



Pr [T r ^T r 



1 + Pr T r = T, 



1 



it follows that E[a(T r )] — rE[Fx r (T r )] > r/2. To prove the sharpness of the bound, consider 

X t =rI(t>U), 

where U ~ U(f3, 1), then T r U and E[a(T r )] = E[rU] = |. To prove it follows immediate by Jenst 
inequality that 

-a(E[T r }) < -E[a(T r )} <-- 



E[T r ] > a" 



One may be interested in obtaining an upper bound for E[a(T r )] or S[T r ] without further assumption. This 
is, however, impossible. The reason is that, without assumptions that control the growth of the process X t 
for any t > T r , the value of a(t) can blow up. Below we introduce a counter-example that demonstrates the 
impossibility that an upper bound can be obtained without further assumption. 
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Example 2.1 Let X t — tY, with Y a non-negative random variable. Then 

a(t) = tE[Y] 

and 

E[T r ] = rE[Y~ 1 ]. 

Suppose Y is exponentially distributed with mean 1, E[Y] = 1 while i?[y _1 ] = oo, so the behavior of a(t) 
is controlled; however, E[a(T r )] = rE[Y]E[Y~ 1 ] = rElY^ 1 ] can be arbitrarily large. Controlling the growth 
of the proces after T r is reached enables us to derive an upper bound as shown in Proposition \ 2. 2\ 

Definition 2.1 A real random variable A is less than a random variable B in the "usual stochastic order" 
if 

Pr(A > x) < Pr(B > x) for all x E (-oo, oo), 
which is denoted A < st B; see \161j and 

Proposition 2.2 Assume that X t is non-negative and continuous, and in addition X t is a time homogeneous 
Markov process and that T^ r — T(u_i\ r > st T r , for k > 2, then 

E[a{T r )] < 2r (6) 

and 

E[T r }<a- 1 (2r). (7) 

Proof 2 First notice that if X t is a Markov process with continuous paths - irreducible state space [0,^4], 
< A < oo. If r < A, then Pr{T r < oo} = 1. Observe that, because of the continuous paths, we have 

T r =T s + (T r - T s ),s < r, 

thus T r is stochastically greater than T s and T r —T 3 . Now 

Pr{T r > x + y} = Pr{T r > x}Pr{T r -x> y\T r > x} 

and 

T r — x\T r > x, X(x) = y ~ T r — T y , 

by the Markov property, and since T r — T z < st T r , for < z < t, Pr{T r — x > y\T r > x} < Pr{T r > y}, so 
that, if we define F(-) = 1 — F(-), 

F Tr (x + y)<F Tr (x)F Tr (y), (8) 

which is the submultiplicative, or new better than used (NBU), property. For details about NBU property, 
see Jlf and J3jj. Since T r is NBU, it has a finite mean. 

Define G(t) := ^ J t °° FT r (x)dx, where fi = E[T r ], the stationary renewal distribution corresponding to 
T r , since 

m = ^ r ^ dx 



V Jt F Tr {t) 

<FtM f°°F T (x-t)dx 
M Jt 

= FtM) 



H = F Tr (t), 



it follows that G < s t Ft t . 



The stationary renewal distribution corresponding to Ft,, has X^ ~ G, and {Xi}i>\ ~ Ft t . It satisfies, 

M*(t) = E[N*(t)] = E[# of renewals in [Q,t]] = -. 

I- 1 
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An ordinary renewal process has X% ~ F with F > st G, since F is NBU. It follows that 

M(t) = E[N(t)} 

= E[# of renewals in [0,t] for ordinary renewal process] 

<E[N*(t)]= t -, 
A* 



id hence 



E[M(T r )] < E\pT l T r ] = 1. 



Under the assumption that T^ r — Tn._\\ r > st T r , then 

N t r < max X s < <N t + l)r, 

0<s<t 

where N t — max{fc : T^ r < t}, i.e. the number of renewals prior to time t. It follows that 

a{t) < rE[N t + 1] 
= (Af ( + l)r 

<(^ + l)r (9) 

and thus 

( E[T, 



E[a{T r )} < + 1 ) r < 2r. 



By plugging in, specifically, t — E[T r ] into we yield 



E[T T ]<a- 1 (^^ + l^r^a- 1 {2r), (10) 



which completes the proof. 

As mentioned, eq. coupled with Eq. (1.1) of [T2"] : 

\a-\ r -)<E[T r ]<a^{2r). (11) 

By assuming that a(t) is concave, the lower bound can be improved to a -1 (|) , which is sharp. If it is further 
assumed that the conditions specified in Proposition 2.2 hold, we obtain the following bounds: 

a- 1 C-)<E[T r ]<a- 1 {2r) 1 (12) 

hence gives the right order of magnitude of the expected value of the first hitting time. In fact, [H] used this 
approach to obtain bounds on -E[T r ] for additive processes including a certain class of stochastic integrals 
extending the works of Burkholder and Gundy; see [5] . 

Readers may compare the results obtained in Propositions 2.1 and 2.2 to Theorem 3 of [T2] whose lower 
bound suggested: E[T r ] > (1 — e)a _1 (er),0 < e < 1, is sharpened by our new lower bound derived in the case 



of concave functions. The results as shown in ( 11 ) and ( 12 ) provide the values of the constants that appear 



in the bounds as shown in |21j . which provides examples of the approach applied to stochastic integrals. 

It can be easily derived from Q and ^ that 

1 < E[a(T r )] - r 

2 - r 



G 



shows stability of E[a(T r )] as r changes and the linearization property of a(t). In addition, this shows the 
linearizing property of a(t) as can be easily seen that by 



sup 

r>0 



E[a(T r )} - 



sup 

!->0 



i?[sup 



<s<T, 



^[sup 



<s<T, 



< 1. 



(13) 



Equation (13) may explain the claim that a(t) can be interpreted as a natural clock for all the processes 
with the same a(t) since through a(t), E[a(T r )] is in a linear relationship with the predehned boundary r. 

The above framework provides a broad foundation for more general applications. A rich array of examples 
are given in [12]. In particular, the above results can be extended easily to X G R d . We will provide the 
following examples for illustration. 

Example 2.2 Suppose that we are interested in the first hitting time of the process X t that hits either a 
lower bound a or an upper bound b, where a < < b. We may define 



fa,b( X ) 



x/a, x < 
x/b, x > 0, 



( see Figure 1 ) and hence 



T a ,b = uif{t > : X t <£ [a, b]} = w£{t > : f a ,b(X t ) > 1} 



a(t) = E 



SUP fa,b{X 8 ) 
0<s<t 



= E 



sup I(X S < 0)a- 1 X s + I(X S > O^X, 

0<s<t 





/ 1 


/ 

/\ 1 
t \* 

f 


1 
1 

1 
1 


□ 


Ta,b 

\ / 


\ 



Figure 1: Illustration of asymmetric bounds for general (not necessarily non- negative) processes. 



More generally, this can be further generalized into cases where the first hitting time is defined as T r = 
mf{t > : f Cl ,c 2 ,....c m (Xt) > r}. Below shown one specical case that can be handled under this framework: 
One can also take for example, (a) f(x) = x + 7 x £ K; (b) f(x, y) = \x — y\,x, y € K d ; (c) f(x, y) = p(x — y) for 
a metric p; (d) f(x,y) = e~\ x ~ y \ , x. y € R d and so on for appropriate applications. Moreover, the boundary 
itself can be not fixed. 

Example 2.3 Suppose that g is a deterministic or stochastic process. Define Y t — X t /gt and T r = vai{t : 
Y t > r}, where X is similarly defined as in Example 3. This corresponds to the hitting time of the process 
X t reaching g t and a(t) — E[swp s<t X s /g s ]. 
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3 Some extensions 



Notice that, under the concave assumption on a(-), we can derive a similar lower bound to the one derived 
in Section 2 that is more readily available. Observe that 



a(t) = E 



sup X s 

,0<s<t 



> sup E[X„] := «(*). 

0<s<t 



It follows that a 1 (r) < k 1 (r). Thus, if the conditions hold for the upper bound, then 

E[T r ] < a" 1 (2r) < K~ X {2r), 
where k(-) may be obtained more easily, compared to a(-). 

Example 3.1 Consider the absolute value of a standard Brownian motion \ W t \. It can be shown that 



(14) 



K (t) = E\W t \ = ^2~tJ^ 

and it is known that -E[sup 0<s<t W s ] = E\Wt\, but a(t) — E[swp Q<s<t \W S \] appears difficult to compute. In 
fact, it equals 

C(t)v^/7T, 

where 1 < c(t) < 2, a multiple depending on t. It should be noted that the conditions for the upper bound to 
hold are easily verifable in this case. Thus 



E[T r ] < K~ X (2r) = 2-nr 2 , 



while it is known that E[T r ] = r 2 



If more information about the process X t is known, we can obtain a relaxed lower bound that can be easily 
expressed in a more manageable form. Suppose X t is a submartingale with right continuous paths, it is well 
known that 

E[X+] V (t) 



Pr< sup X s > r } < 

[0<s<t J T T 

Therefore, it follows that Pr{T r < t} — Pr{sup 0<;s<t X s > r} < r)(t)/r, which leads to 

E[ v (T r )] 



> J Pr{T r < t}dF Tr (t) > X - 



E[v(T r )} > 



(15) 



which is a stonger result compared with Proposition 2.1. 

Example 3.2 Consider W t , where Wt is a standard Brownian motion. Denote T r the first passage time of 
\Wt\ to r, (so for the first passage time to r 2 ). Applying the result of (15), we get 

E[T r ] > y. 

Note that the actual value of E[T r ] is r 2 in this case. Here, rj(t) = E[X 2 ] = t. 



The usefulness of (15) can be demonstrated in the following example in which a closed form expression of 
E[T r ] is difficult to obtain. 

Example 3.3 Consider a submartingale 

X t = \Wf-t\, 

where Wt is a standard Brownian motion. It can be shown that 



V (t)=E\W 2 ^t\=tE\W 2 -l\ 
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Then it follows that 

r 



so that 



2 <E[ V (T r )] = ^-E[T r ] 



T / 7T6 

E[T r ] > - J— «0.5116r. 



In this exmple, the mean ofT r , the first passage time of \Wf — t\ to r, is difficult to compute. 

Before ending this section, we would like to point out that the finiteness of the mean of T r is important 
because this issue can cause problems in other applications. A good illustration is demonstrated in the 
following example. 

Example 3.4 Consider 

W t + = vaax{0,B f }, 

where, again, B t is a standard Brownian motion. The first passage time of to r > coincides with the 
first passage time of Wt to r which has an infinite mean. W t + is not a Markov process. Hence, results of 
Proposition 2 do not apply to case. T r is finite with probability 1, but this is not NBU; see definition 

Since W s +t ~ W s + {W s +t — Wt) — W s + Wt, where Wt is an independent Brownian motion, 

W£ r+t ~ {W Tr + W t ) + = {r + W t ) + <r + W t + . 

When W^~ hits r, +r hits 2r but W^ +t may not have hit 2r. Thus T 2r > s t T r + T r and T 2r — T r > T r . 
More generally, T^ k+1 y — T)- r > st T r . Thus, following our previous argument, we have 

rN t < max W s < r(N t + 1), 

0<s<t 

but in this case, {N t } is a renewal process with an infinite mean interarrival time. There is no stationary 
distribution G (of T r ) and ([9| does not hold. In this case, E[a(T r )] = E[*J 2/7P\/!Zr] = oo. This shows that 
upper bounding E[a(T r )] is challenging to deal with in general. 



4 The rate of growth of the maximum of Bessel processes 

This section is dedicated to the discussion of the rate of growth of the maximum of Bessel processes which 
can be obtained via the inequalities obtained in Section 2 (and 3). We are going to consider a case in which 
the hitting time of the radius/surface area/volume of the largest multi-dimensional Brownian motions that 
hits a predefined boundary. 

Example 4.1 Let pi = (j>\ ,p\ ,p\ ), i — 1, ...,m, be a set of m spheres that perhaps represent some 
identified tumours in a human body, i.e. a three-dimensional space. At time t = 0, we start a three- 
dimension Brownian motion with coordinates processes centered at each one of these points. For each i a 
sphere of radius r*j (t) is given, where the radius equals the distance of the location of the three dimensional 
Brownian motion at time t to its starting point We are interested in getting qualitative information on 
how long it will take before the radius (volume) of at least one of the spheres exceeds a fixed level (say r, 
r > 0) as the size of m varies. 

Let X = {Xi(t) = (XP(t),XP(t),X^(t))} i=h 2,..., m where, for each i, X;(*) = p s ; + B((t) corre- 
sponds to the location of the three-dimensional Brownian motion Bj(t) started at point Pi and B(t) = 

i=i,...,m- The radius of the largest sphere is given by \\B(t)\\ = 
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Figure 2: Illustration of Examples 5 and 6. 



Observe that for each j, Yj(t) — Y^_^(Bj (t)) 2 is a submartingale, since 

E[Y 3 (t)\F s ] = j^E [{Bf {t)f\F s ] > (s)f = Yj(s). 



i=l 

It follows that Y t :— maxi<j< m Y (t) is a submartingale because 

E[Y t \T s ] = E[max{Yi(t), . . . ,Y d {t)\T s }] 

>max{E[Y j {t)\F s ],j = l,...,d} 
> max(Yj(s), j = 1, . . . ,m) = Y s . 



Let T r .d be the first passage time of \\B(t)\\ = sup J<m yXii=i(^i''W) 2 t° r ' equivalently the first passage 
time ofY t = \\B t \\ 2 to r 2 . Then 



Pr{T r!d < t} = Pr |supF s > r 2 | 



< E W 1 p 



max xl, 

Ki<m 



Therefore, we can write 

1 < E [maxi< t < d xi,i] 

2~ r 2 ■ - l- - ,»j - 2E [maxKKd X3 



- ^2 ^[^dj => #Pr,<iJ > - - r ^^r. (16) 



can fee shown that E [maxx<i<d xl »] — ^Si=i ' 1 a7lc ^ hence (16) can &e rewritten 



(IS 



E[T r A > f . (17) 

which follows from the Corollary on Page 266 of f7]/. In fact, we can also approximate the values of 
-E[maxi<i<d X3 %] y i a simulations. The corresponding values for various d are tabulated in Table [7J' 

In general, for {Xi(t),t > 0, }i=i,..., m be iid as {X(t),t > 0} a submartingale. Let be the first passage 
time to r for {Xi(t),t > 0}, then mini<i< m T^ is the first passage time to r for maxi<j< TO Xi(t). Let 
{X(t),t > 0} be independent of {Xi(t),t > 0}i = i ... ;m with T r , its first passage time to r. Note that 

J Pr{T r < t}dF inin{T W i ... iT w } (*) = Pr {T r < min^, . . . , T^}} > (d + l)" 1 . 
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d 


1 2 3 4 5 10 


E 


y/maxi<;< d xi,i 




1.599 1.979 2.173 2.324 2.413 2.720 


d 


15 20 30 40 50 100 


E 






2.875 2.987 3.132 3.237 3.307 3.527 



Table 1: Approximation of E ^/maxi<i<dx|^ via 10,000 simulations 



Since 

p r{ T r <* } <ra = ^, 

r r 

so 

(d+^r 1 < J Pr{rr<0^ min{ T 1 W,..., m i„ { rM}(*) 

^KminlT^...,^ })] ^ 
r 

As a result, ^^(minjT^, . . . , T^})} > ^ For ||A t ||, we use || -Bj || and r 2 instead. Recall that 



E 



ELil^i W) 2 =3*. This gives 



min Ti 

!< 3 <d 



> 



3(d+ 1)' 



This result is not as good as what we can obtain in Example 3.4 in which E[T r ,d] > „„J — — 7. But there 
we used the result that %| is IFR (increasing failure rate); in other examples, we might know the type of 

(r) 

distribution that T> follows. 

The following example studies the upper bound of the Bcsscl process discussed in Example 3.4. 

Example 4.2 (Example 5 of a non-Markov process for which the upper bound holds) Consider the Bessel 
process studied in Example 5 again. \\B S \\ is known to be strongly Markov but maxi^x^ ||-B^|| may not be. 
Note that 

Bt+s = B s + (B t+S — B s ) =d B s + B t , 
where {B t ,t > 0} is distributed as {B t ,t > 0} and B is independent of B. It follows that 

\\Bt+ s \\ =d \\B S + B t \\ < st \\B S \\ + \\B t \\, 

and hence 

\\B t +s\\ - \\B S \\ < st \\B t \\, t>0, 

independently of T s , the history accumlated up to time s. Thus, independently of T s , the first passage time 
of 1 1 .Bt+ S 1 1 — ||B S || to r is stochastically larger than T r . 

Now consier \\{B^} t \\, j = 1, . . . ,d. Given T s , which denotes the history of all j processes up to time 

s, the conditional distribution of each ||Bt+ s — -Bg^|| is stochastically larger than each of the \\B^\\ for all 

t > 0. It follows that the first passage time of Z t := maxi^jx^ ||Bj || from kr to (k + l)r is stochastically 
larger than the minimum of d random variables, each distributed as T r . Hence T? k+1 -. r — T£ r > st T* for all 
k > 1 . It further follows that, T kr is stochastically larger than the convolution of n iid random variables, 

each distributed as T* , where T* = mini^xd (= T r ^). 
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Letting N t = m&x{k : T* r < t}, N(t) = max{k : J^Li T ri < t }> 

we can write 
1 ' + 1 ) r 



1 



EN t < E[N(t)] < 

\E[1 r ,d\ 

since T* — T r> d is NBU. Next, denote Yd = d in&xi<i/<<i{X3 u }, we have 



V~tE[^/¥ d ] = E 



max IIS; 

i<j<d" 



< r 



t 

E[T*] 



It follows that yjE \T r ]E[^JYd\ < 2r [by letting t = E\T Tid \ and \jE\T r ] < 2r(E\/Y d )]. As a result, we have 

4r 2 



E[T r ] < 



[EVYd? 



(18) 



The two-sided bound in this case is thus 



< E[T r>d ] < 



4r 2 



or equivalently, 



A[E^T d ] 2 ~ {EVYd}*' 
a- 1 (r/2)<E[T r ^ d ]<a- 1 {2r), 



(19) 



which is the same as what we obtained in eq. (12 1. Note that, in this example, the underlying process 
Z t = maxi<j-<d H-Bp^ll is not Markov. However, each \\B^\\ is Markov and so T r (r) is NBU. It follows that 
T* = minKjXd T$?' is NBU and thus has a finite mean. 

Finally, we would like to emphasize that since the constants involved in the bounds derived are independent 
of the size of d, the inequalities obtained can be used to derive quantitative comparisons on the expected first 
passage times for processes with different values of d. That is, if we include the dependence on d for different 
values, say di and d 2j an d take ratios, we have 

mTP _ E[su Vi < di IBM)*] 
E[sup,< d2 \By\p\ 

which gives us information on the relative growth rate between the maxima of Bessel processes. 

All the examples shown in this section involve Brownian motion, below shown is an example that demonstrate 
how the bounds can be applied to other types of random variables whose distribution is not Gaussian. 

Example 4.3 Let X±, X2, ■ ■ ■ , X n be non-negative, possibly dependent random variables with Pr{Xj 
A" > = for all i ^ j. Let T r denote the rth smallest amongst X\, . . . ,X n , Fi be the marginal CDF of 
X t and N(t) = {ffX, < t}. Then, 



E[a(T r 



(*)=.E[iV(i)] 

i=l 

n 



E 



> 



Thus, N(T r ) = r and 



E[a(T r )]=E[N(T r )} > -, 



illustrating the decoupling aspect of this inequality. Again, N(t) = {ffXi < t}, where {Xi}i = i^^ n are inde- 
pendent copies of {-2Q}i=i,...,« It should be noted that this lower bound is obtained without the knowledge of 
the dependence structure of X 's. 
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The above setting can be used to model a pool of debtors whose survival times (time until which they 
become default) follow some distribution with non-increasing density function, say exponential distribution 
or a subset of Weibull distribution. In this case, T r can be interpreted as the time when r (out of n) debtors 
have gone bankrupt, which can be an important time stamp that triggers the termination of payment to a 
lower tranch of collateralized debt obligation (CDO). By assuming stationarity, we can treat the data as two 
sets of independent copies and use the historical data to estimate the current set of individuals (random 
variables). Unlike the use of copula to model the event time, the results presented previously can provide 
bounds for the event time without knowing the dependence structure of the debtors. 



5 Conclusion 

In this paper, we derive bounds for the expectation of the stopping time of arbitrary stochastic processes. 
The approach we use (see [7] and |12j ) involves the concept of boundary crossing of a non-random function 
to that of a random function. In the situations where the moment of the maximal process is available, the 
results shown can be helpful for the estimation of E[T r ]. In particular, for non-negative, continuous and 
time homogenuous Markovian processes, with the assumption that T^ r — T\k-i)r ^st T r for k > 2, we show 



that the order of magnitude is the same up to constants. The result of (11) suggests that it is appropriate 



to view the process through a(t), which serves as a clock for all processes with the same a(t), as reflected in 



(13 1. The lower bound derived can be applied to arbitray measurable processes and it is particularly useful 



in the study of renewal processes. 
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